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ABSTRACT 

System identification techniques are used to 
represent the dynamic behavior of the Space Shut- 
tle Main Engine. The transfer function matrices 
of the linearized models of both the closed loop 
and the open loop system are obtained by using 
the recursive maximum likelihood method. 

NOMENCLATURE 

A , B , C state space matrices 

a denominator polynomial of the transfer 

function matrix 

b numerator polynomial of the transfer 

function matrix 

D signal duration, sec 

e equation error 

f frequency, rad/sec 

H,Hjj system transfer function matrix and its 
i , j element 

MR oxygen to fuel mixture ratio of the 

main combustion chamber 

n discrete time integer 

N number of samples 

P pressure, psi 

S speed, rpm 

T temperature, °R 

t time, sec 

t c clock time, sec 

u.x.y input, state and output vectors 

z Z transform variable 


♦Summer Faculty Fellow in 1988; presently 
at Mechanical Engineering Dept., Florida 
Atlantic University, Boca Ratan, Florida 33432. 


Greek Symbols 


Su deviation of input from nominal 


Sx 

deviations of state from nominal 

Sy 

deviations of output from nominal 

0 

valve actuator output rotary motion 
dimensionless 

Subscr i 

£tS 

CCV 

coolant control valve 

C 

chamber 

FPOV 

fuel preburner oxidizer valve 

HPFT 

high pressure fuel turbine 

HPOT 

high pressure oxidizer turbine 

max 

maximum 

mi n 

minimum 

MFV 

main fuel valve 

MOV 

main oxidizer valve 

OPOV 

oxidizer preburner oxidizer valve 

R 

request 


INTRODUCTION 


Accurate representation of the dynamic 
behavior of the Space Shuttle Main Engine (SSME) 
is required for a variety of diagnostic and con- 
trol design and evaluation purposes. A complete, 
nonlinear dynamic simulation has been developed 
(Rocketdyne Division of Rockwell International 
Corporation, 1981). However, its size and com- 
plexity make it very difficult to use for design 
purposes. This paper describes the identifica- 
tion of linearized dynamic models of the SSME 
from the nonlinear dynamic simulation. The iden- 
tified linearized models are valid in limited 
response regions about the operating point cor- 
responding to the 100 percent power level. They 
are suitable for the purpose of control system 
design and development. 


1 


Initially a brief description of the SSME 
is given. This is followed by a description of 
the identification process and the model used. 
Finally, some results obtained from the identi- 
fied models are compared with the results 
obtained from the nonlinear simulation for the 
same input. 

THE SPACE SHUTTLE MAIN ENGINE 

Different aspects of the SSME as well as 
its principle of operation are described in the 
literature (Landauer, 1988; Klatt et al., 1982). 

A brief description of the main engine is given 
below by following them. 

The space shuttle orbiter main propulsion 
system is composed of three space shuttle main 
engines. The engines use liquid oxygen and liq- 
uid hydrogen propellants that are carried in an 
external tank attached to the orbiter. These 
elements and the two strap-on solid propellant 
rocket motors are assembled in the vertical posi- 
tion. (Solid rocket boosters and the SSMEs pro- 
vide the thrust for the first portion of the 
flight until the solid propellant motors burn out 
and are separated.) The SSME ' s continue to burn 
until the vehicle is near the desired orbital 
velocity. Following SSME shutdown, the external 
tank is separated and final orbit insertion is 
accomplished by the orbital maneuvering system. 
After completing orbital operations, the orbiter 
reenters the atmosphere and glides to a landing. 

The engines operate for a total of 480 sec 
from launch. Each produces a sea level thrust 
of 375 000 lb and a vacuum thrust of 470 000 lb. 
To understand the overall flow of fuel and oxi- 
dizer to produce the thrust, a schematic diagram 
of the propellant flows and the control valves 
i s shown in Fig. 1 . 

Pressurized fuel, provided by the fuel tank 
flows through the low pressure fuel pump and the 
high pressure fuel pump, is fed to the regenera- 
tive cooling and the preburners. A pressurized 
oxidizer tank provides the oxidizer. It flows 
through the low pressure oxidizer pump, to the 
high pressure oxidizer pump where the output 
flow splits as shown in Fig. 1. 

The two high pressure turbines are driven 
by a fuel turbine preburner and an oxidizer tur- 
bine preburner, each of which produces hot gas. 
The low pressure turbines are driven from the 
high pressure pump flows. The fuel from the 
high pressure fuel pump (HPFP) goes through the 
main fuel valve (MFV). After the MFV, the flow 
divides into fixed nozzle cooling flow, main 
chamber cooling flow, and coolant control valve 
(CCV) flow. The relative proportions of these 
three flows is determined by the CCV. Heat is 
absorbed from the combustion chamber nozzle and 
the fuel flows to the preburners where combustion 
and pressure is controlled by the fuel preburner 
oxidizer valve (FPOV) and the oxidizer preburner 
oxidizer valve (OPOV) . In the fixed nozzle cool- 
ing, the coolant control valve determines the 
percentage of the fuel used for cooling the com- 
bustion chamber by establishing a bypass flow 
direct from inlet to discharge. 

The high-pressure fuel rich hot gas from 
the preburners goes through the two hi gh-pressure 
turbines and returns through the hot gas mani- 
folds to the fuel injector (FI) of the combustion 
chamber. Main chamber cooling fuel flow powers 
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the lew pressure fuel turbine (LPFT) and returns 
directly to the fuel injector of the main combus- 
tion chamber. 

Oxidizer from the high pressure oxidizer 
pump (HPOP) splits into three paths. Flow from 
the first path drives the low-pressure oxidizer 
turbine (LPOT). The output flow from the LPOT 
merges with the output of the low pressure oxi- 
dizer pump (LPOP). The second flow goes through 
the main oxidizer valve (MOV) to the oxidizer 
injector (01) in the combustion chamber. The 
third flow enters the preburner oxidizer pump 
(P80P) and is divided between the preburners as 
determined by the FPOV and the OPOV. This com- 
pletes the flow paths for turbine drive and 
control . 

THE CONTROL SYSTEM OF THE SSME 

The existing control system of the SSME uti- 
lizes five valves FPOV, OPOV, MFV, MOV, and CCV 
to control the mixture ratio and the main chamber 
pressure which is correlated with the thrust. 

Open and closed loop control of these valves are 
used to accomplish the SSME mission. A typical 
SSME mission is shown in Fig. 2. 

All the valves are initially scheduled to 
follow the start sequence as indicated in Fig. 3. 
During the main stage, the thrust and the mixture 
ratio are controlled by adjusting the OPOV and 
FPOV valves. This adjustment is based on thrust 
and mixture ratio request, through a closed loop 
proportional plus integral feedback control. A 
conceptual block diagram of the control system 
during the main stage of operation is shown in 
Fig. 4 . 

For the purpose of model identification the 
following inputs and outputs are used. The 
closed loop system inputs are the mixture ratio 
request, MRr, and the chamber pressure request, 
Pcr, which is proportional to the thrust. The 
open loop system inputs are the angular displace- 
ments of the valve actuators for each of the 
five valves, ©OPOV* ©FP0V> ©CCV. ©M0V> ©MFV • 

The outputs are the chamber pressure (P(0, HPFT 
pressure (P^p ft), HPOT pressure (PhpoT^* HPFT 
temperature (Thpft>, HPOT temperature (Thpq|), 
mixture ratio (MR), HPFT speed (Shpft^* and HPOT 
speed (Shpot^- In this study all parameters are 
normalized and expressed as the fractional devia- 
tions from the nominal operating condition. 

SYSTEM IDENTIFICATION 

Identification may be defined as the deter- 
mination of the mathematical model of a process 
among a class of mathematical models based on 
the measurements of input and output data of the 
process. Hence, the steps in identification con- 
sist of: determining a priori knowledge, selec- 

tion of a driving signal that contains enough 
spectral components to excite all inodes of the 
system, selection of a model, and parameter esti- 
mation. These steps are followed in the identi- 
fication of SSME dynamics and are outlined below. 

A Priori Knowledge 

The responses of both the closed and the 
open loop system for test input signals are exam- 
ined to obtain preliminary information about the 
types of nonlinearities and the bandwidth of 
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these systems. The responses of the closed loop 
system and some of the selected outputs of the 
open loop system as well as the test input sig- 
nals are shown in Figs. 5 to 7. 

The closed loop response, as shown in 
Fig. 5, indicates existence of both the stiction 
and backlash type of nonlinearities and a small 
system bandwidth of approximately 5 Hz. The 
backlash nonlinearity is more obvious from the 
response of the open loop system, since the out- 
puts cannot return to their initial values when 
the step input is removed, as shown in Fig. 6. 
Figure 7 indicates that CCV, MOV, and MFV valve 
actuator angular displacements are essentially 
decoupled from the outputs. Therefore at this 
nominal operating condition only the actuator 
angular displacements of the valves FPOV and OPOV 
are used as the inputs of the open loop system. 

Selection of a Driving Signal 

In this study, a pseudo random binary 
sequence (PRBS) is selected as the input per- 
turbation signal of the system because of its 
convenience and its suitability for similar 
applications (Cottington and Pease, 1978, 

Sudhakar et al., 1988). A PRBS is a wideband, 
long-duration signal which possesses one of two 
values. The signal switches from one value to 
another at time intervals determined by a clock. 

A PRBS is easily generated using shift registers 
or a computer program, and can be made small 
enough amplitude so as not to significantly 
alter the operating condition of the engine and 
yet still yield useful information for identifi- 
cation. The amplitude of the PRBS signal used 
in this work is 2 percent of the nominal operat- 
ing point magnitude of the chosen input. 

The clock time, t c , is chosen to be in the 
range of system time constants as follows. The 
minimum and the maximum observed frequencies, 
f m i n and f max , are estimated by examining the 
response of both the closed loop and the open 
loop system from Figs. 5 and 6. Then the clock 
time, t c , and the number of samples needed are 
determined by using the estimated values of the 
minimum and the maximum frequencies using the 
equations below: 


D = J — (2) 

min 


where D is the signal duration and N is the 
number of samples. For the identification of 
the open loop system, a clock time of 0.04 sec 
and 511 samples are used. This corresponds to a 
maximum frequency of 78.5 rad/sec (12.5 Hz), a 
minimum frequency of 0.31 rad/sec (0.05 Hz) and 
a signal duration of 20.44 sec. A step input is 
used as the driving signal for the closed loop 
system since its response is very slow, as shown 
in Fig. 5. 


Selection of a Mathematical Model 

As mentioned earlier, Figs. 5 and 6 indicate 
the existence of significant backlash and stic- 
tion nonlinearities associated with the valve 
dynamics. It is assumed that these two nonlin- 
earities can be isolated from the rest of the 
system as shown in Fig. 8. Here, the open loop 
system includes actuator dynamics, the nonlinear 
element which contains valve stiction and valve 
linkage backlash, and the engine system which 
includes all other nonlinearities and the engine 
dynamics. To simplify the analysis, the backlash 
and stiction were removed from the simulation. 

All subsequent results are obtained with these 
two nonl i neari ties removed. Now the open loop 
system, without the backlash and the stiction, 
was linearized and identified. The backlash and 
the stiction nonlinearities can later be added 
to the identified system. In practice values for 
backlash and stiction would be determined from 
manufacturer ' s specifications or bench testing 
of the individual components. 

The open loop dynamics of the SSME can be 
described by the nonlinear equations: 

x(t) = f[x(t),u(t>] (4) 

y(t) = g[x( t) ] (5) 

where x, u, and y are the state, the control, 
and the output vectors. Linearizing these equa- 
tions about a nominal operating condition and 
discretizing yields: 

Sx(n + 1) = A6x(n) + BSu(n) (6) 

Sy(n) = CSx(n) (7) 

where Sx, Su, and Sy are the deviations of the 
state, the input and the output vectors, respec- 
tively, about the nominal operating condition. 
Taking the z-transform of Eqs. (6) and (7), the 
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system model relating the input vector 

Su 

and the 

output vector Sy can be obtained as: 



Sy(z) = H(z)Su(z) 

(8) 

where 




H. J < 2 > = C. (zl - A) -1 Bj 

(9) 


Here, C \ is the i row of the C matrix and 
Bj is the j th column of the B matrix. A 
parameter estimation algorithm can be used to 
obtain the parameters of each of the elements of 
the transfer function matrix, Hij(z). 

Parameter Estimation Algorithm 

The recursive maximum likelihood ( RML) algo- 
rithm (Eykhoff, 1981) is used to estimate the 
parameter values. It computes the coefficients 
of the transfer function between an input and an 
output. That is, if 6ui<n) and Sy](n) are the 
input and output of a linear time invariant sys- 
tem, H] i (z) , then the RML algorithm determines 
H] i ( z) from the data Su](n) and Sy](n) for some 
range of time n. H]j(z) is found by minimizing 
the total energy of the so called equation error: 
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e(z> = a(z) 6 y 1 ( 2 ) - b(z) 6 u^ (z) 


( 10 ) 


H(z) 


b(z) 

a(z) 


( 11 ) 


Generally, the RML algorithm needs fewer 
data samples or parameter estimation than other 
approaches and can account for the presence of 
certain types of noise. It can be modified to 
track slowly time varying parameters and can be 
extended to the analysis of multi-input multi- 
output systems (Eykhoff, 1981, Ljung and 
Soderstorm, 1 985) . 

"he RML method estimates the parameters of 
a moael when its order is known. However, the 
determination of the order of a model is still an 
open problem (Ljung and Soderstorm, 1985). The 
method employed here is based on the "Parsimony 
principle." In this method parameters are esti- 
mated for models of increasing orders. The esti- 
mate of the final error norm is computed for 
each, and the order for which the decrease In 
error is no longer appreciable is adopted (Box 
and Jenkins, 1970). Thus, from the error plot, 
which plots error versus order of the model, the 
optimum order is chosen as the one corresponding 
to the "knee' of the error plot. A typical 
error plot is shown in Fig. 9. 

The output of the nonlinear simulation for 
PRBS'input perturbations is used to determine the 
transfer function matrix of the linearized model 
by utilizing the RML algorithm. Each transfer 
function matrix element is determined by utiliz- 
ing the response of the system to small perturba- 
tions about a nominal operating point as follows. 
The element of the input vector is per- 

turbed and applied to the system with other ele- 
ments of the input vector set at their nominal 
values to obtain the perturbed outputs. By sub- 
tracting the nominal values from the perturbed 
values, the differential input Suj(n), and the 
corresponding outputs Syi<n), for i =1,2, 3, 
... are obtained. This procedure is repeated 
for al 1 the i nputs . 

The RML algorithm outlined before is avail- 
able in a software package called Ctrl-C (Systems 
Control Technology, 1986), and is utilized to 
estimate the parameters of the model. Since the 
nonlinear simulation is used as the unknown sys- 
tem, the system is assumed to be noiseless. 

Under this condition the RML method simplifies 
to the recursive least squares method. 

RESULTS 

The procedure outlined earlier is used to 
identify both the closed loop and the open loop 
system. The responses of the identified model 
are compared with the responses obtained from 
the nonlinear simulation for the same input and 
observed to be in good agreement as shown in 
Figs. 10 and 11. The parameters of the elements 
of the transfer function matrix are given in 
Table 1. 

The responses of the open loop system are 
also compared with the responses obtained from 
the nonlinear simulation. A typical result is 
shown in Fig. 12. In general, it is observed 
that the gains for positive and negative pertur- 
bation signals are significantly different. In 
order to compensate for this phenomenon, a sys- 
tem gain with different values for positive and 


negative perturbations is added to the linearized 
model. In fact it is this ratio of the magni- 
tudes of the positive and negative values which 
must be identified. This ratio is equal to the 
ratio of the magnitudes of the positive and nega- 
tive steady state values of the outputs obtained 
from the nonlinear simulation. One model that 
includes this variable system gain can be 
expressed mathematically as: 

Sy(z) = H(z) 6 u*(z) (12) 

aSu(k) Su(k) > 0 

Su*(k) = (13) 

(2 - a) 6 u(k) ou(k) < 0 

The identification procedure is carried out 

with this modified PRBS and a piecewise linear 
model is obtained. Comparison of the responses 
of this model with the responses of the nonlin- 
ear simulation showed excellent agreement as 
shown in Figs. 13 to 16. The parameters of the 
coefficients of the transfer function matrix as 
well as the values of a are given in Table 2. 
The proposed model block diagram structure is 
given in Fig. 17. 


CONCLUSION 

System identification techniques are used 
to represent the dynamic behavior of the SSME. 

The comparison of the responses of the nonlinear 
simulation with the responses of the identified 
model indicates very good agreement. The identi- 
fied model can be used for control design 
purposes . 

The identified model does not include valve 
linkage backlash and valve stiction nonlineari- 
ties. These nonlinearities should be added to 
the identified model and the validity of the 
model should be checked by comparing it with the 
full nonlinear simulation. 

The identified model is valid for a limited 
response region about the 100 percent power level 
operating condition. This study will be extended 
so as to obtain models at other operating condi- 
tions. Then these point models can be linked to 
cover the full range of operation of the SSME. 
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TABLE 1. - IDENTIFIED MOOEL PARAMETERS 
FOR THE CLOSED LOOP SYSTEM 

+ b-z ~ 2 + b^z ” 2 + b.z ” 4 

( 7 ) _ ! - i Z 

ij 2 ”, -1 -2 -3 

J 1 + a 1 z + a 2 z + a^z + a 4 z 


Input 

P CR 

MRr 

Output 

p C 

MR 

Pc 

MR 

Parameter 







d l 

-1 .048 


3 


) 

-1 .025 

d 2 

.285 





.325 

d 3 

.102 





-.223 

a 4 

-.134 





0 

b] 

-.008 





.017 

b 2 

.272 





.336 

b 3 

.053 





-.276 

b 4 

-.113 





0 


TABLE 2. - IDENTIFIED MODEL PARAMETER FOR THE OPEN LOOP SYSTEM 

b,z -1 + b ? z “ 2 + b-.z ” 3 + b.z “ 4 

H..(z> = j—* “3— 1 — -7 

J 1 + a^ z + a 2 z“ z ♦ a^z J + a 4 z 4 


(a) Input = ©OPOV 


Output 

Pc 

P HPFT 

p HP0T 

t HPFT 

t hpot 

MR 

Shpft 

Shpot 

a 

0.88 

0.88 

0.89 

0.85 

0.94 

0.9 

0.87 

0.89 

d l 

-.676 

-.708 

-.787 

-.821 

-.457 

-.526 

-.998 

-.501 

d 2 

.0373 

.0602 

.0985 

.16 

.0151 

.04 

.252 

-.023 

d 3 

0 

0 

.0183 

.0151 

.0448 

0 

0 

0 

bl 

.0517 

.1093 

.5358 

-.117 

2.911 

.0555 

.0092 

.2164 

b 2 

.4947 

.4627 

.0086 

.914 

-.429 

.4975 

.1416 

.4698 

b 3 

0 

0 

.0925 

-.533 

-.004 

0 

0 

0 


(b) Input = ©FPOV 


Output 

Pc 

p HPFT 

p HP0T 

t HPFT 

t hpot 

MR 

S HPFT 

Shpot 

a 

N/A 

0.8217 

N/A 

0.8707 

0.9681 

0.9882 

0.896 

N/A 

41 


D 

-1 .043 


9 

-.604 

-.272 

-.237 

-.393 

0 

d 2 



-.062 



.075 

-.046 

-.007 

-.028 



d 3 



.211 



0 

0 

0 

0 



bi 



.3419 



1.502 

-1 .001 

-.115 

.0885 



b 2 



-.181 



-1 .223 

-1.312 

-.858 

.2742 



b 3 



-.08 



0 

0 

0 

0 
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SSME PROPELLANT FLOW SCHEMATIC 
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FIGURE 1. - PROPELLANT FLOW AND CONTROL VALVES (ROCKWELL, 1983). 



FIGURE 2. - TYPICAL 104% SSME MISSION (ROCKWELL INTER- 
NATIONAL CORP., 1988). 
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V START ENABLE COMMAND (ENGINE READY SELF MONITOR) 
V PRE-CHILL BLEED VALVES CLOSED 


V ENGINE START COMMAND 
JvMAIN FUEL VALVE STARTS RAMP TO FULL OPEN 
J V FUEL PREBURNER OXIDIZER VALVE STARTS TO OPEN 

] V OXIDIZER PREBURNER OXIDIZER VALVE STARTS TO OPEN 

I 

| V MAIN OXIDIZER VALVE STARTS OPEN AT 60% 

[ y SCHEDULE PREBURNER 

^ VALVES FOR PRIMING 

| V SCHEDULE CHAMBER COOLANT VALVE 

I THRUST LOOP CLOSED 

j V THRUST TO 100% POWER LEVEL 

I V RAMP MAIN OXIDIZER VALVE FULL OPEN 

, MIXTURE RATIO LOOP CLOSED V 



“ 2-10 1 2 3 4 5 

TIME FROM ENGINE START COMMAND, SEC 

FIGURE 3. - SSME START SEQUENCE (ROCKWELL INTERNATIONAL CORP., 1988). 
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COMMAND CHAMBER 



ENGINE 

OXIDIZER 

FLOWRATE 




MIXTURE RATIO = 


CONSTANT X CHAMBER PRESSURE 


FUEL FLOWRATE 


FIGURE 4. - SSME PERFORMANCE CONTROL (ROCKWELL INTERNATIONAL CORP., 1988). 


7 




NORMALIZED OUTPUT 





9 







NORMALIZED OUTPUT NORMALIZED OUTPUT 



(a) MULT I PULSE RESPONSE. 


ORIGINAL 
OF POOR 


p AGE i S 
Q' W JTY 



FIGURE 11. - COMPARISON OF THE RESPONSES OF THE 
LINEAR MODEL WITH THE NONLINEAR SIMULATION 
(CLOSED LOOP). INPUT = MR R ; OUTPUT = MR. 
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o (a, CHAMBER PRESSURE OUTPUT. 

LU 



(b) HIGH PRESSURE OXIDIZER TURBINE PRESSURE OUTPUT. 


FIGURE 12. - COMPARISON OF THE RESPONSES OF THE 
LINEAR MODEL WITH THE NONLINEAR SIMULATION 
(OPEN LOOP). INPUT = 0 opov ; OUTPUT = S HppT . 


FIGURE 13. - COMPARISON OF THE RESPONSES OF THE IDEN- 
TIFIED MODEL WITH THE NONLINEAR SIMULATION (OPEN 
LOOP), (a) INPUT = 0 opov ; OUTPUT = P c< (b) INPUT 
0 opov ; OUTPUT = P Hp0T . 
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NORMALIZED OUTPUT 
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FIGURE 16. - COMPARISON OF THE RESPONSES OF THE 
IDENTIFIED MODEL WITH THE NONLINEAR SIMULATION 
(OPEN LOOP), (a) INPUT = 0 Fp0V ; OUTPUT = MR. 
(b) INPUT - 0 Fpov ; OUTPUT = S HpFT . 
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FIGURE 17. - PROPOSED SSME SYSTEM DYNAMICS MODEL AT 100% RATED POWER. 
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